Geometric phase of an atom in a magnetic storage ring 
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A magnetically trapped atom experiences an adiabatic geometric (Berry's) phase due to changing 
field direction. We investigate theoretically such an Aharonov-Bohm-like geometric phase for atoms 
adiabatically moving inside a storage ring as demonstrated in several recent experiments. Our result 
shows that this phase shift is easily observable in a closed loop interference experiment, and thus 
the shift has to be accounted for in the proposed inertial sensing applications. The spread in phase 
shift due to the atom transverse distribution is quantified through numerical simulations. 

PACS numbers: 03.65.Vf, 39.20.+q, 03.75.-b, 39.25.+k 



INTRODUCTION 



A recent highlight in atom optics ll[ is the realiza- 
tion of atomic storage rings 0, IH, H, U 0] . With well- 
designed magnetic (B-) field configurations, cold atoms 
including Bose-Einstein condensates can be confined in a 
ring shaped region with a diameter of the order of a mm 
0, a cm 0, or 10 cm 0,[{| 0. Such a trapping geometry 
provides exciting prospects for delivering atoms into high 
Q optical cavities 0] and for inertial sensing applications 



Magnetic ring shape traps are formed by inhomoge- 
neous magnetic fields. Accompanied by the trapped cen- 
ter of mass motion, the atomic spin adiabatically fol- 
lows its local B-field direction, giving rise to a geomet- 
ric (Berry's) phase 0, |, [T|. This induces a phase 
shift in storage ring based interference, reminiscent of 
the Aharonov-Bohm effect from the magnetic gauge po- 
tential in electron interference. Even when confined in 
transverse motional ground state, a distribution of this 
geometric phase arises corresponding to atomic guiding 
trajectories at different transverse radii, potentially caus- 
ing degradation to the interference contrast. Therefore a 
detailed investigation is essential for applications of ring 
based interferometry to high precision measurements. 

The geometric phase of an atom in the widely used 
Ioffe-Pritchard trap was first studied in Ref. , where 
the associated gauge potential (Berry's connection) was 
also discussed for a general magnetic trap. The aim of 
this paper is a detailed investigation of this geometric 
phase for an atom in a magnetic storage ring. We find 
that it is simply related to the cosine value of the an- 
gle between the local B-field and the symmetry axis of 
the ring, a result easily understood in terms of the solid 
angle sustained by the adiabatic changing B-field along 
the central azimuthal closed trajectory of a storage ring. 
We have performed extensive numerical studies of the 
geometric phase as a function of the various components 
forming the trap field and elucidated the condition when 
the fluctuation of this geometric phase is small. 

This paper is organized as follows: in Sec. II, we briefly 
review the atomic storage ring system and discuss the 



associated geometric phase. Based on the assumption of 
a simple classical center of mass motion of the atom, the 
geometric phase and the spatial fluctuation are calculated 
in Sec. Ill for various types of magnetic storage rings. 
In Sec. IV, we provide an full quantum treatment of 
this geometric phase in a hypothetical storage ring based 
interference experiment. Finally we conclude in Sec. V. 



II. GEOMETRIC PHASE IN A STORAGE RING 

As in recent experiments H, 0, [E 0j a magnetic 
storage ring can be formed by several axially-displaced 
concentric current carrying coils. Near the origin, the 
trap field takes the form B = B p e p + B z e z [HI, HH with 
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up to polynomial terms second order in the cylindrical 
coordinate (p,<j),z). The B\ terms clearly resemble the 
quadruple B-field from a pair of anti-Helmholtz coils, 
while the Bq and Bi terms take the familiar expansions 
for a Ioffe-Pritchard trap [HI]. It is easy to verify that 
when BqB2 > 0, the trap field vanishes at the point 
Poe p + z e z [11( of 
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-B x jB 2 
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To avoid spin flipping (Majorana transition) near this 
zero B-ficld region, an azimuthal bias field 0] can be 
added, or a time averaged potential can be constructed 
by rapidly oscillating the trap field B [ll[. In this 
paper we mainly focus on the more complicated case of a 
time averaged trap (TORT) from the recent experiment 
0. When -Bo, 1,2 are oscillating functions of time with 
a large frequency u>, the atomic motion will be governed 
by a time averaged potential cx (oj/2tt) \B(r,t)\dt, 
provided adiabaticity is maintained. This requires u> to 
be much less than the atomic Larmor frequency in the 
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trap field such that its spin remains aligned. Also u> has 
to be much larger than the time averaged trap frequency 
in order to make the trapping effective. 

For simplicity, we will suppress the (p, z) degrees of 
freedom and consider atomic motion parameterized by 
<f> along the azimuthal direction of the storage ring. 
Then the state of the atomic mass center would be de- 
scribed by the function <f>(t). The interaction between 
an atom and the presumably weak trap field B[<f>, t] is 
governed by the linear Zeeman effect 

H[fat] = n B g F P • B[<f>,t], (3) 

where pb is the Bohr magneton, and F is the atomic 
hyperfine spin. The Lande factor for the spin-1 manifold 
of a 87 Rb atom is cjf = —1/2. For simplicity we take H — 
1. The trapped atomic spin is adiabatically aligned with 
respect to the local B-field, i.e., the atomic spin remains 
in the low field seeking internal state |(— l)[</>, £])b of the 
instantaneous eigenstate of H with eigenvalue e_i[f] = 
—Pb9f\B\- It is related to the eigenstate |(— 1)) 2 of F z 
(of eigenvalue —1) by a simple rotation 

|(-1)[<M)b - exp[-iP • e,/3(i)]|(-l))», (4) 

where (3(t) is the angle between B[<j>, t] and the z-axis. It 
is independent of <fi due to the cylindrical symmetry. 

In this section and the following section III, we treat 
the atomic center of mass motion is classically, by its 
one dimensional azimuthal angle <f>(t) along the circum- 
ference of the ring (at fixed p = p c and z = z c ). The 
quantum treatment of the atomic center of mass motion 
will be given in section IV. According to the quantum 
adiabatic theorem [3], an atom initially described by 
state |W(0)) = |(-1)[0(O),O]) S becomes 

|*(t)>= |(-l)[^),t]) B e-^«e-^( t ) ! (5) 

at time t. 7d = J Q * e-x[t']dt' is the dynamical phase, and 
7 s (i) is the geometric phase, i.e., the Berry's phase [|[ 
that is expressed as 

iS) = -i j\{{-i)[mA\^\{-i)W')A)Bdt' 

= f cos[(3(t')}n{t')dt' , (6) 
Jo 

with £l(t) = d(j)/dt the atomic angular velocity 
along the storage ring. To obtain the above for- 
mula we used d/dt = d t + {d<p/dt)d^, as well as 
the equalities l)[^,t]|9^|(— l)[(/),t\)B = «cos/3 and 
s((— 1)[0, i]|9t|(— 1)[0, <])b = that come directly from 
our definition of |(-T)[</>, t]) B in Eq. 

The periodic time dependence of B[<f> } t] gives 
rise to the Fourier expansion cos[/3(t)] = cos/3o + 
Sn>o cos[nu>t + (p n ] with cos/3o denoting the time av- 
eraged value of cos[/3(t)]. On substituting into Eq. ©, 
we find 7 9 = cos /3 - <£(0)] + J2n>o C ™ I cos[nut' + 



ip n ]Q,(t')dt' . The integral C n J cos[nujt' + (p n ]n(t')dt' 
takes the maximum value 2C n tt / (nu>) for a slowly varying 
function f2(£), which becomes negligible when f2 is much 
smaller than u>. In this case, we find the geometric phase 
7 9 can be approximated as 7 9 (t) ~ cos /3q [<p(t) — </>(0)], 
which is independent of £l(t) and completely determined 
by the trajectory of the atom in the storage ring. For a 
closed path, we arrive at the geometric phase 

7c = 2tt(7Cos/3o, (7) 

where q is the familiar integer winding number of the 
atomic trajectory. More specifically, we find 

. w [ 2 * /uJ B z {r c ,t) 
cospo = — / — dt, (8) 

within our approximation for various trap field configu- 
rations. In the above, f c — p c e P + z c e z denotes the center 
position of the time averaged storage ring trap. 

III. RESULTS 

In the previous section, the general expression for the 
geometric phase ([7]) is obtained for an atom in a storage 
ring. This section is devoted to a study of the size and 
distribution of this phase for various types of magnetic 
storage ring. For a static storage ring trap formed by 
supplementing an azimuthal bias B-field from a current 
carrying wire along the symmetric axis as in Ref. [6| , the 
Berry's phase vanishes because /3o = ir/2. For a TORT 
trap, the dependence of cos /3q on the trap parameters 
-80,1,2 of Eq. {1} can be investigated. When B2 is a con- 
stant, cos/3o is closely related to the oscillation amplitude 
of Bi . For the following simple time dependence 

B a = B 2 [L 2 + n 2 sin(wt)] , B 1 = B 2 l cos(ojt) , (9) 

parameterized by length scales L, n, and I, we find that 
the value of cos/3o is mainly determined by l/L. In Fig. 
[IJa), the l/L dependence of cos/3o is shown for different 
values of n/L. We see that cos/?o behaves as an increas- 
ing function of l/L. For n = 0, it turns into an approx- 
imate linear dependence cos/?o = 0.31/ L. In Fig. [1Kb), 
the storage ring trap center p c is shown for n = 0. It 
behaves like a quadratically decreasing function of l/L. 

During oscillation, the zero field point of Eq. ([2]) traces 
out a closed trajectory in the cross section of the storage 
ring. In the proposal of Ref. [ll[ and the experiment 0], 
this closed trajectory forms a loop (a torus in 3D) sur- 
rounding the center (p c , z c ) along the storage ring. For 
a constant B 2 , the oscillation amplitude of B\ is pro- 
portional to the radius of this "death torus" @ in the 
^-direction. If the radius becomes sufficiently small, e.g., 
as in Fig. 3(b) of Ref. [llj, the value of cos/?o in the 
trap center also becomes very small. The locus of the 
zero point for our example is plotted in Fig. HJa). The 
trajectory that forms the boundary of the loop (torus) 



FIG. 1: (Color online) (a) The dependence of cos/3o on l/L for 
n — (solid line), n = 3L/4 (dashed line), and n = L (dash- 
dot line), (b) The dependence of pc/L on l/L for n — 0. 




p/L p/L 



FIG. 2: (Color online) (a) The locus of the zero field point for 
the trap's B-field. (b) The equal B-field contours of the time 
averaged magnetic field. The B-field is in units of Bo. Both 
figures are for n = and I — L. 

has now collapsed into an open curve on one side of the 
trap center, although the time averaged trap potential is 
still well established as illustrated in Fig. E^b). We have 
also calculated the Fourier coefficients C n of cos @(t). In 
Fig. 3(a), the two largest coefficients Ci and C4 are plot- 
ted as functions of l/L. We note that their maximum 
values remain much smaller than unity. 

In the above calculation, the classical center of mass 
motion of the trapped atom along the storage ring is at 
a fixed r c . In reality, the atomic center of mass posi- 
tion (ft, z) may deviate from (p c , z c ) even for a closed 
trajectory. This uncertainty of atomic position in the 
cross sectional plane can be due to either transverse 
atomic thermal motion or the quantum nature of the 
transverse motional state. In the following calculation 
for the fluctuations, we consider the simple case of a 
transverse distribution reflected for instance by its mo- 
tional ground state. This then leads to a genuine de- 
pendence of the geometric phase on the values of trans- 
verse coordinates p and z. This causes fluctuations of 
the phase shift, which can lower interference contrast or 
even destroy interference fringes from the atomic wave 
packet in the <j> direction as we discuss in next section. 
Within the semi-classical framework we adopt, the fluc- 
tuations can be approximately simulated by a proba- 
bilistic sampling of the values for p and z as reflected 
by their probability distributions in the transverse mo- 
tional ground state. To arrive at a high contrast in- 
terference pattern, the spread of this phase fluctuation 
should be kept small. We have calculated the varia- 
tions of this fluctuation with respect to different param- 



FIG. 3: (Color online) (a) The l/L dependence of C 2 (C 4 ) in 
blue (red) squares, (b) The fluctuation of cos /3o over the re- 
gion 5 = 0.015L (blue empty square), S = 0.005L (red circle) 
and 5 — 0.001L (black plus sign). 

eters of a TORT trap. In Fig. [3jb), the fluctuation 
defined according to / = ([cos Po(r) — cos (3o(r c )] 2 ) is 
plotted as a function of l/L for 6 = 0.001L, 5 = 0.005L 
and 5 — 0.015L. Here cos @o(r) is the time averaged 
value of B Z /\B\ at the position r = pe p + ze z . The 
average (•) is over a properly weighted, assumed here 
for simplicity as a flat spatial distribution in the region 
R : {p € [p c -S,pc + 8],ze [z c -5,z c + 6}}. f is found to 
be a decreasing function of l/L. Together with Fig. [lja), 
we conclude that when l/L is large, the value of cos/3o is 
also large with small fluctuations. 

Finally we illustrate a practical example for B2 = 7800 
Gauss-cm -2 and L = I = n = 0.1 cm. The bias field 
Bo(t) is 78[1 + sin(wi)] Gauss and the first order B- 
field gradient B\(t) becomes 780sin(w£) Gauss-cm -1 . In 
this case the trap frequencies in the p- and z- directions 
are 345Hz and 672Hz, respectively. The time averaged 
Berry's phase is 2ir cos /?o ~ 0.2-7T. This trap can be re- 
alized with two pairs of Hclmholtz coils (a and b) cen- 
tered at (z = ±j4 a ,b, p — 0) and a pair of anti-Hclmholtz 
coils c centered at (z = ±A C , p = 0) with A a — 0.1cm, 
Ab = 0.25cm, and A c = 0.5cm, respectively. The radii 
R a .b,c of the three pairs of coils are taken as R a = 0.3cm, 
Rb = 0.5cm, and R c = 0.6cm with currents in the two 
Helmholtz pairs being i a = 289A and ib = — 550A, and 
in the anti-Hclmholtz coils are i c — 335A. 



IV. THE AHARONOV-BOHM LIKE PHASE 
SHIFT FROM A QUANTUM TREATMENT 

The geometric phase as discussed above is directly 
measurable, e.g., in the Sagnac interferometer [7]. For 
this purpose, we perform a calculation capable of a 
quantum treatment of atomic motion. Trapped atoms 
are assumed to be in the guiding regime so that their 
transverse degrees of freedom, p and z, are frozen in 
the ground state. Writing the atomic wave function as 
t)\(—l)(cj>, t)) b with t) for its spatial motion, we 
find that the Schrodinger equation i(d/dt)^> = -ffad^ is 
governed by a time averaged adiabatic Hamiltonian [IBJ 

#ad = ^ [Pcb-A(cf>)/ Pc } 2 + V (<£). (10) 
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V{4>) denotes the time averaged trap potential 
that more generally can contain additional opti- 
cal or the gravitational potentials. A{<j>) = 

~i(tu/2n)f^ /uJ B <(-l)(0,t)|0*|(-l)(&t)> fl dt is the 
time averaged gauge potential. For the simple case con- 
sidered previously, A equals cos (3q and is independent of 
<f>. In this section, without loss of generality, we assume 
A(<p) to be a slowly varying function of <j>. The operator 
P4> = —ipc 1 d/d t p denotes the atomic canonical momen- 
tum along the direction. In the presence of the gauge 
potential A(<p), the atomic kinetic angular velocity 
becomes v$ = p$/ (mp c ) - A{<j>)/ {mp 2 c ). 

We note that the Hamiltonian in Eq. (flU)) only ap- 
plies when the gauge potential A do not depend strongly 
on p and Z . clS otherwise their dependence in the term 
A((j>, p, z)p<f, would lead to inseparable correlations be- 
tween the dependencies on cf>, p, and z of the atomic 
wave function. If this does occur, the reduced atomic 
quantum state for the cj> direction becomes a mixed state 
rather than a pure one, resulting in reduced coherence 
property. This constitutes a quantum explanation for 
the reduced contrast in the interference pattern due to 
fluctuations of the geometric phase. 

Obviously this Hamiltonian (TlU|) resembles the mo- 
tion of an electron around a current carrying solenoid, 
where an Aharonov-Bohm (A-B) phase [161 ] is induced 
by the corresponding vector potential, with a value pro- 
portional to its integral along a closed path surrounding 
the solenoid. In our case, we expect essentially the same 
result with an Aharonov-Bohm-like phase, or the Berry's 
phase 7c as defined in Eq. J7]), given by the line inte- 
gral of the induced gauge potential J A(<fi)d(j>. In the A-B 
effect of an electron, both the electron source and the de- 
tector (screen) are far away from the solenoid, where the 
vector gauge potential vanishes. The electron achieves 
a steady asymptotic distribution long before and after 
(t — > ±oo ) the encounter with the solenoid. Therefore its 
motion can be treated via stationary scattering theory 
0, E3] or based on a path integral approach with an in- 
cident plane wave [l8| . For our case of trapped atoms in 
a storage ring, the gauge potential A((p) never vanishes at 
any instant during the evolution. Thus our model must 
be treated dynamically as described by a time dependent 
Schrodinger equation. 

We first consider the motion of a single wave packet 
assuming that A(<fi) varies slowly and is approximated as 
a constant in the region, where the wave function is non- 
negligible. The expectation value of the atomic position 
cj), the canonical momentum p^, and the velocity of 
the wave packet at time t will be denoted as <j>t, pt, and 
v t . At t = the wave packet y(</>, 0) can be written as 



¥>(</>, 0) = V(</>,0)e iA( * o)( *-* o) . 



(11) 



Both fa and vq are determined completely through 

V>W>,o), 



\\f_e 



ii(0>|> 




FIG. 4: (Color online) A proposed setup of interferes in a 
storage ring. An initial superposition of two wave packets 
^ ± (</,y A(0W with counter propagating velocities ±uo inter- 
ferences at time T when the packets overlap near the other 



end 



"o 



mp c 



0)cV#, 0)#, (12) 



independent of the gauge function A. 

The Hamiltonian H a d (flU)) can be obtained from the 
one without gauge potential A(<f)) via a transformation, 
i.e., H ad = UHoU' 1 with 



Ho = + V(<f>), U = exp[z f 4 ' A{cj>')d^'] 



(13) 



provided the wave packet is limited to and moving in a 
finite region over the circumference. Therefore, (p((p,t) is 
written as 

^(<M) « ^(0 > t) e "te>W-&> e i '& A( * W , (14) 

with ip((f),t) — e~ lHot ip((f>,0), where we have used the 
relation e~ lHt — Ue~" lH ° t U~ 1 and replaced A{4>) with 
A{4>q) \A{<j>t)] m the region where tp(cj), 0) [tp((j),t)] is non- 
negligible. The center position (velocity) (pt (vt) of each 
wave packet and the variation of the profile ip are also 
determined by Hq independent of A(<p). 

To detect the geometric phase, we consider an inter- 
ference set up as shown in Fig. [4] with the initial atomic 
state being a superposition of two counter propagating 
wave packets ip±e lA ^ > ^ centered at <f> = 



$(0,0) 



1 

V2 



iM<^)e^ (o)0 ± ?M0, 0)e lA(0 H • (15) 



The two wave packets overlap again at time T when 
<^ +) = 0^ + 2?r. Making use of Eq. (14]), we find 

$(^ T ) = ^[V, + (0,T)e 4A ^ +) )^-^ +) )e 4 ^ +> ^W')<i0' 
v2 

+V-(^- 27r J T)e iA ^ +> )(*-4 +) ) e «/o^ +) ~ i " r A(*')<i*'] j (16) 



o = / </#W>,0)| 2 #, 



with tp±(<f),T) = e~ lHo1 ip±(4>, 0). A substitution of <j> -> 
cf>— 2ir to the counter propagating wave packet defines the 
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complete wave function over the same azimuthal region 
(0, 2tt). An interference pattern oc cos[r] + ((fi) — ?7-(</>)+7c] 
then reveals the presence of the phase shift 7c, where 77+ 
(r;_ ) is defined as the phase of ifi+ (</>, T) [ifi- (<fi — 2tt, T)] . 
Since the evolution of ifi± is governed by Hq , the density 
distribution n±((j>) for each component and the phase 
difference Tj+((fi) — V—ifi) are independent of A(<fi). The 
induced gauge potential causes a shift, simply the geo- 
metric phase shift 7c, to the interference pattern. 

In many situations, the center of mass motion for the 
two interfering paths of an atom can be described clas- 
sically. In this case the interference pattern can be pre- 
dicted with more detail. We show in the appendix that 
Eq. (fT4|) . which describes the motion of a single wave 
packet, can be generalized to the following 

<p(<f>,t) « e l ^[^>*^o,O] e4 [p t -po][0-0 t ] e -»[<? t -<?o]p^(0 ]O ), 

(17) 

with fa and pt denoting the expectation values for the 
position and its canonical momentum that satisfy the 
Hamiltonian equation governed by -ff a d, i-e., the center 
velocity v t of the wave packet satisfies 

v t =Pt- A(4> t ) = dfa/dt. (18) 

As we have pointed out before, the classical trajectory 
(fit and the speed Vt are independent of the gauge po- 
tential A((fi) because the dynamics is one dimensional. 
The global phase S c \ is defined as the action function in 
classical mechanics: 

S* c i(^,t;0 o ,O) = / L[4> v ,v v ]dt', (19) 
Jo 

where L\<fi t , v t ] = v\ /2+A{<fi t )v t — V{<fi t ) is the Lagrangian 
function in classical mechanics. With this semiclassical 
description for the wave packet inside the storage ring, we 
find that the relative phase [r]+((fi) — r}-{4>) + 7c] can be 
approximated more transparently as [mp c (v^ —v^, )(fi-\- 
£ + 7c] with a gauge-independent constant 

£ = r ?0 +(0)-7 70 -(0)+ / [(p<4 +) ) 2 /2-y($->)]df 
Jo 

- f T [(pA^Y/i-v^Vm. (20) 

Jo 

The discussions outlined in this section also allow for 
easy calculations of more general Aharonov-Bohm like 
phase shifts due to different gauge effects in storage ring 
based matter wave interference experiments. For ex- 
ample, the matter wave Sagnac effect, which can be 
considered as an analog of the Aharonov-Bohm effect 
is easily deduced to correspond to a phase shift 
2^rr£lp 1 c cos 9 in a storage ring set up, with f2 being the 
angular velocity for rotation and 9 being the angle be- 
tween the axis of rotation and the z direction. The deduc- 
tion simply involves replacing the adiabatic gauge pro- 
tectional A(<fi,t) with mf^cos^, which is proportional 
to the </>-component of the " gauge field" mVt x r for the 
Sagnac effect |19j ]. 



V. CONCLUSION 

In this paper, we have presented the theoretical anal- 
ysis of the geometric phase for a neutral atom inside a 
magnetic storage ring. For the interesting example of 
a time averaged potential (TORT), we have shown that 
the geometric phase is proportional to the averaged co- 
sine value of the angle between the local B-ficld and the 
z axis. When the oscillation amplitude of the quadruple 
component B\ is large, a significant and sharply peaked 
geometric phase is realized on completing a single pass 
along the storage ring. Of course the net effect is cumula- 
tive with respect to the number of turns an atom makes. 
We hope our result will shine new light on the proposed 
inertial sensing experiments based on trapped atoms in 
storage rings. 
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APPENDIX A: THE ONE DIMENSIONAL 
SEMI-CLASSICAL MOTION OF A SINGLE 
WAVE PACKET 

During the course of this work, we obtained some inter- 
esting results on the semi-classical one-dimensional mo- 
tion for an atomic wave packet. These are summarized 
here because to our knowledge, they do no seem to have 
been stated or known anywhere. 

Essentially the following material provides a proof of 
Eq. (|17p which describes the one dimensional semi- 
classical motion of a single wave packet tp(<fi,t). Our 
calculation is based on the following two assumptions: 

1. During the whole evolution, the wave packet is as- 
sumed to be sufficiently narrow that the potential 
V(<fi) can be approximated by a linear function in cfi 
and the gauge potential A^ assumed a constant in 
the region where the wave packet is non-negligible; 

2. The evolution time is so short that the effect of 
wave packet diffusion can be omitted. 

With the first assumption, it is easy to prove [2(| that 
the center position cfi t , the center canonical momentum 
(fit, and the center velocity Vt satisfy the classical Hamil- 
tonian equations 

d ^ - - A/1 \ 

-TT = v t =p t - A((fit), 
at 

f - < A1 > 

As we have shown in Eq. (HH), at t = 0, the wave 
packet can be written as (p((fi,0) = ifi{cfi, Q) e iA ^ a ^-^ . 
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Furthermore, we can express the profile ij){4>, 0) as 

ip{<p,0) = I dpf(p,0)e- ip ^-^e- ao< -' t> -^ \ (A2) 



where the function f(p, 0) describes the momentum dis- 
tribution at t = 0. 

At time t, the atomic wave function ip(<p, t) is expressed 
as in Eq. ([H]). determined by = e~ iHot ^p(<j), 0). 

To calculate the wave function i/j(4>,t), we first replace 
the trap potential V{4>) with a time dependent linear 
potential Vo(t) + Vi(t)(f> with 



V (t) = V{4> t )- 

dV{4>) 



dV(4>) 



<t>=<t>t 



Vi(t) = 



<t>=<t>t 



The Hamiltonian equation (jAlj) satisfied by the center 
position and velocity can now be casted in the form of a 
Newtonian equation 



d 2 4> f /dt 2 = -Vx(t), 
dcj)t/dt = v u 



(A3) 



which can be solved as 

4>t — 4>o + v t 



t ft! 

dh / dt 2 Vi(t 2 ), 



t't 



JO 

dhVi(h). 



(A4) 



According to the first assumption, the wave function 
ijj((j),t) — e~ lHo ijj(({>,0) satisfies the Schoredinger equa- 
tion 

idtVOM) = HesiWifat), (A5) 
approximately, with the effective Hamiltonian 



H eS (t) 



V Q (t) + Vx{t)<i>. 



(A6) 



Because the effective Hamiltonian H e g is a quadratic 
function, the Schrodingcr Eq. (|A5I) can be solved with 
the Wei-Norman algebra method [2l|. After direct cal- 
culation, we find 

i>(<P,t)=K(t)iP(<j>,0), (A7) 

with the evolution operator given by 

U{t) = e-^He-^^e-^^e-^, (A8) 

parameterized by 



o(t) = t/2, 
b(t) = /* dhhVxih) 
c(t) - fidhViih), 

d{t) = /*dt 1 Vb(ti) + Jo^i/o 1 dt2Vi(ti)t2Vi(t 2 ) 
In the end, we find 



(A9) 



^(0,t) = e iSo(t) / dp/(p,i)e- ipW -* t) e- lCtW -^ ) , (A10) 



with /(p, t) = /(p, 0)e *p 2 '/ 2 and the global phase factor 



S (t) 



L [<f>t',v t >]dt', 



(All) 



with L [4> t ,v t ] = Vf/2— V a (t) — Vi(t)4>t- In the above 
derivation, we have used Eq. (| A4|) and the relationship 
f*dt 1 {t-t 1 )V 1 {t 1 ) = f*dh f* 1 dt 2 {t-t 1 )V 1 {t 2 ). 

If t is so small that the momentum uncertainty Ap of 
the wave function "0(^j 0) remains small, we can omit the 
wave packet diffusion and assume /(p, t) « /(p, 0), which 
is our second assumption above. With this approxima- 
tion and substituting Eq. (|A10[) into Eq. (fT4)) . we obtain 
Eq. (jTTJ) . Of course, this same result can also be deduced 
with the method in Sec. 3 of Ref. 221. 
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